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A legend tells |l|] that once Loschmidt asked 
Boltzmann on what happens to his statistical 
theory if one inverts the velocities of all parti- 
cles, so that, due to the reversibility of New- 
ton's equations, they return from the equilib- 
rium to a nonequilibrium initial state. Boltzmann 
only replied "then go and invert them". This 
problem of the relationship between the micro- 
scopic and macroscopic descriptions of the phys- 
ical world and time-reversibility has been hotly 
debated from the XlXi/i century up to nowadays 
|2|^]. At present, no modern computer is able to 
perform Boltzmann's demand for a macroscopic 
number of particles. In addition, dynamical chaos 
1 10 13] implies exponential growth of any impre- 
cision in the inversion that leads to practical ir- 
reversibility. Here we show that a quantum com- 
puter [14 -17 1 composed of a few tens of qubits, 
and operating even with moderate precision, can 
perform Boltzmann's demand for a macroscopic 
number of classical particles. Thus, even in the 
regime of dynamical chaos, a realistic quantum 
computer allows to rebuild a specific initial dis- 
tribution from a macroscopic state given by ther- 
modynamic laws. 

To study the relations between microscopic determinis- 
tic classical dynamics, macroscopic thermodynamic laws 
and quantum computation, we choose a simple area- 
preserving map: 



dw(y,t)/dt = D/2 d 2 w(y,t)/d 2 y, 



(2) 



y = y + x (mod L) , x — x + y (mod 1) 



(1) 



Here the first equation can be seen as a kick which 
changes the momentum y of a particle, while the second 
one corresponds to a free phase rotation in the interval 
—0.5 < x < 0.5; bars mark the new values of the vari- 
ables. The map dynamics takes place on a torus of integer 
length L in the y direction. For L — 1 this map reduces 
to the well-known Arnold cat map Jl(|, which describes 
a fully chaotic dynamics with positive Kolmogorov-Sinai 
entropy h w 0.96. As a result, the dynamics is character- 
ized by exponential divergence of nearby trajectories, so 
that any small error e (for example round-off error) grows 
exponentially with time, and reversibility of a trajectory 
is lost after ig w |lne|//i map iterations. For e ~ 10~ 8 
comparable to ordinary precision of the Pentium III, this 
time scale is rather short (tg 20). For L 1 chaos 
leads to the diffusive spreading of particles in momentum, 
which is well described by the Fokker-Planck equation: 



where the diffusion coefficient D «< x 2 >= 1/12. Thus 
after a time t > an initial distribution of particles 
in (Q) evolves towards a Gaussian statistical distribution 
w(y, t) = w g (y, t) = exp(-(y - y ) 2 /(2Dt))/ V2Dnt with 
< y 2 >= Dt + y 2 , where y =< y > at t = 0. On a fi- 
nite torus this diffusive process relaxes to a homogeneous 
distribution in y after a time to ~ L 2 / D. 

For the case L — 1 it was shown that a quantum com- 
puter can simulate a discretized version of this map with 
exponential efficiency jl8| . Here we show that for L> la 
similar quantum algorithm enables to simulate the evolu- 
tion of a macroscopic number of classical particles which 
is governed by the thermodynamic diffusion law. To per- 
form this evolution on a lattice of size LN 2 (with TV = 2 TLq 
and L = 2 riq '~ riq ) this algorithm uses three quantum reg- 
isters. The first one with n q qubits holds the values of 
the coordinate x (a;, = — 0.5 + i/N,i = 0,...,N— 1), 
the second one with n q i qubits holds the y coordinates 
(ijj = -L/2+j/N,j = 0,...,LN-1) and the last one with 
n q i — 1 qubits is used as a workspace. The first two regis- 
ters describe the discretized classical phase space with L 
cells and ./V 2 points per cell. In this way, the initial po- 
sitions of Nd ~ N 2 particles can be represented by one 



quantum state ^ 



. a ij\ x i > IVj > |0 >i where a itj = 
or 1/y/Nd- The quantum algorithm is based on modu- 
lar additio ns p erformed in a way similar to the one de- 
scribed in fl9[ , through Toffoli and controlled- not gates 
(CNOT). It requires 10n g + 6n g ' — 17 gate operations per 
map iteration, in contrast to 0(2 2riq ) operations for the 
classical algorithm. The time inversion is also realized by 



8n 9 + 4n, 



ql — 13 gate operations which effectively change y 



into — y half-way between kicks |1S| [E0[. In this way, the 
quantum computer acts in a way similar to Maxwell's de- 
mon |2l],|22| who reverses the velocity of each individual 
particle. 

A perfect quantum computer simulates exactly the 
map (^) , but a realistic physical system always has some 
imperfections which can destroy time-reversibility. For 
a classical computer, e. g. Pentium III iterating map 
(|l|), round-off errors of amplitude e destroy the time- 
reversibility of the map dynamics after iterations. 
This fact is illustrated on Fig.l where it is assumed that 
the demon inverts the velocities of all trajectories after 
t r = 35 iterations with a precision e. After that, the 
macroscopic distribution starts to return back but after 
tg ~ \ hie\/h iterations the errors become too large and 
the diffusion process restarts again. In contrast to that, 
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on a quantum computer round-off errors can be decreased 
enormously since the size of the registers grows only lin- 
early with n q and the computation is exponentially effi- 
cient. However, a quantum computer has its own natural 
errors which can be viewed as imprecisions of amplitude e 
in the gate operations. The comparison of the two types 
of errors natural for classical and quantum computers is 
displayed on Fig.l. It shows that the quantum compu- 
tation with precision e = 0.01 in each gate at each map 
iteration is able to reverse effectively the diffusion pro- 
cess up to the initial state. That is in striking contrast 
with the irreversibility of the classical computation with 
a round-off error of amplitude e = 10 -8 made only once 
at t = t r when the demon acts. In this way, the quantum 
computer succeeds to reverse the thermodynamic diffu- 
sive process with enormous number of particles. Indeed, 
at the moment of inversion t r , the distribution of parti- 
cles is a Gaussian of width a = \J2Dt r in agreement with 
the solution of (||), as is shown in Fig. 2. 
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FIG. 1. Diffusive growth of the second moment < y > of 
the distribution w(y, t) generated by the map (|l|) with L = 8, 
simulated on a classical (Pentium III) and quantum ("Quan- 
tium I") computers. At t = t r — 35 one inverts all velocities. 
For Pentium III inversion is done with precision e = 10 -4 
(red line) and e = 10 -8 (green line); 10 6 orbits are simu- 
lated, initially distributed inside the demon image (see Fig. 
3). For Quantium I, the computation is done with 26 qubits 
(n q = 7, n q i = 10) (blue line) ; each quantum gate operates 
with imperfections of amplitude e = 0.01 (unitary rotation on 
a random angle of this amplitude). The black straight line 
shows the theoretical macroscopic diffusion with D — 1/12. 

Fig. 3 shows explicitly the distribution in phase space 
at different moments of time. The initial distribution 
mimics a demon, which at t = t r is transformed to a sta- 
tistical homogeneous distribution in the x direction, with 
a smooth variation in y described by (Q). The quantum 
computer operating with 1% accuracy is able to recover 
the initial image with good precision, whereas the classi- 



cal computer with round-off errors 10 -8 completely fails 
to reproduce it. The striking difference between the two 
final distributions at t = t2 r generated by the two com- 
puters can be easily detected from a polynomial number 
of measurements. 
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FIG. 2. Distribution of particles in y for map (gj) simulated 
by Pentium III (green lines) and Quantium I (blue symbols) 
for the case of Fig.l, at t = 20 (dashed lines and open circles) 
and t — t r = 35 (full lines and filled circles). Black lines show 
the theoretical solution of the Fokker-Planck equation (|2j). 



The previous results are supported by the data for the 
fidelity f(t) defined as the projection of the quantum 
state in presence of gate imperfections on the exact state 
without imperfections. For / = 1 both states coincide, 
whereas for / -C 1 both distributions are completely dif- 
ferent. The results on Fig. 4 show that f(t) smoothly 
decreases with number of iterations t even if classical 
dynamics is exponentially unstable. The probability of 
transition from the exact state to other states induced 
by imperfections can be estimated as of the order of e 2 . 
Hence, since imperfections in each gate are assumed to 
be uncorrelated, f(t) should drop by n q e 2 at each map 
iteration (for n q ~ n q i). This determines a time scale 

t f « C/(n q e 2 ) (3) 

on which the fidelity of quantum computation for the 
algorithm is reasonable (/(*/) = 0.5)), even in absence 
of error correction. This scaling is in agreement with 
the data in Fig. 4 (see also jjl|) which give the numerical 
factor C w 0.5. This is in sharp contrast with classi- 
cal errors for which computation of trajectories remains 
correct only up to a time scale tg w |lne|//i. It is in- 
teresting to note that the situation is similar to the time 
evolution of a physical system in the regime of quantum 
chaos, which is stable against small quantum errors even 
though the underlying classical dynamics is chaotic p3| . 
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FIG. 3. Evolution of a demon image through map (|l|). Left 
column shows the simulation on Pentium III, right column 
on Quantium I. Top: initial distribution in the central cell 
(-0.5 < x,y < 0.5). Middle: distribution at t = t r = 35 
in the whole phase space (—0.5 < x < 0.5, —4 < y < 4). 
Bottom: distribution at t — 2t r = 70 in the two central cells 
(—0.5 < x < 0.5, — 1 < y < 1). The time-inversion is made at 
t r = 35, with accuracy e = 10~ 8 for Pentium III (error is done 
only at t r ), and with accuracy e = 0.01 for Quantium I (error 
is done at each gate operation). Color marks the density of 
particles/probability, from blue (minimal) to red (maximal 
value). Here as in Fig.l n q — 7,n q i = 10, with in total 26 
qubits used for Quantium I; for Pentium III, 10 6 orbits are 
simulated. 
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FIG. 4. Fidelity / as a function of e n q t for the quantum 
simulation of the map (|l|) on Quantium I. Here L = 8 and 
n q — 7, (black line and symbols), n q = 6 (blue symbols), 
n q = 5 (red symbols), n q = 4 (green symbols)and e = 10~ 2 



(diamonds), e = 3.10 
up). Black line shows n q = 7, e = 10 
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The relation (||) implies that a quantum computer op- 
erating with realistic accuracy can invert velocities of all 
particles at a given moment of time, so that a specific 
initial state is reliably reproduced from a statistical dis- 
tribution described by a diffusive process. Such a sim- 
ulation for a macroscopic number of particles can 
be performed with few tens of qubits. For example, for 
N d = 6.022 x 10 23 (Avogadro's number) the simulation 
with L = 8 requires only 125 qubits. Moreover, accord- 
ing to (JsJ) , the accurate simulation of such an enormously 
large number of particles remains reliable (f(t) > 0.5) up 
to a time t w 150 with a modest gate accuracy e = 0.01. 
Such a computation is far out of reach of any modern 
supercomputer, and clearly shows the power of quan- 
tum computers. It also opens interesting perspectives 
for cryptography since an initial image can be coded in 
a thermodynamic distribution with very large entropy 
and then reliably recovered. Thus quantum computers 
open new possibilities for the investigation of the rela- 
tions between microscopic deterministic dynamics and 
macroscopic thermodynamic laws. 
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